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Abstract 

An analysis of the Independent Ramjet Stream (IRS) 
cycle is presented. The IRS cycle is a variation of the con- 
ventional ejector-ramjet, and is used at low speed in a rocket- 
based combined-cycle (RBCC) propulsion system. In this 
new cycle, complete mixing between the rocket and ramjet 
streams is not required, and a single rocket chamber can be 
used without a long mixing duct. Furthermore, this concept 
allows flexibility in controlling the thermal choke process. 
The resulting propulsion system is intended to be simpler, 
more robust, and lighter than an ejector-ramjet. The perfor- 
mance characteristics of the IRS cycle are analyzed for a 
new single-stage-to-orbit (SSTO) launch vehicle concept, 
known as “Trail blazer.” The study is based on a quasi-one- 
dimensional model of the rocket and air streams at speeds 
ranging from lift-off to Mach 3. The numerical formulation 
is described in detail. A performance comparison between 
the IRS and ejector- ramjet cycles is also presented. 

Introduction 

The NASA Glenn Research Center is currently develop- 
ing a reusable, single-stage-to-orbit (SSTO) launch vehicle 
known as “Trailblazer 1 ” that utilizes a rocket-based com- 
bined-cycle (RBCC) propulsion system. This vehicle, shown 
in Fig. I, will operate in four modes from lift-off to orbit; 1) 
air-augmented rocket, 2) ramjet, 3) supersonic combustion 
ramjet (scramjet), and 4) all-rocket. The Trailblazer is an 
axisymmetric vehicle having three RBCC propulsion pods 
equally-spaced at \ 20 ° intervals. This configuration insures 
that the total thrust vector is always aligned with the vehi- 
cle's axis, and allows the forebody boundary layer to be 
diverted between the pods. An axisymmetric vehicle was 
chosen in part, for its high structural and volumetric effi- 
ciency, and for its low drag, at the cost of lower inlet pre- 
compression. The entire cross-sectional area of the vehicle is 
available for nozzle expansion. A detailed description of the 
Trailblazer vehicle is given in Ref. 1 . 


A propulsion system using the ramjet and scramjet cycles 
must have another means of acceleration from static condi- 
tions to low supersonic speeds, at which point the ramjet 
cycle can generate sufficient thrust for further acceleration. 
RBCC systems use rocket motors to accomplish this, and are 
characterized by a high degree of integration between the 
rocket and ramjet cycles. This integration can provide both 
thermodynamic and structural advantages over a system in 
which the rocket and ramjet engines are separate. Some form 
of RBCC propulsion is generally considered to be an appro- 
priate choice for air-breathing, reusable, SSTO launch 
vehicles 2 . 

In this study, we focus on the low speed mode which typ- 
ically covers the speed range from lift-off to a maximum 
around Mach 3. Two general classes of RBCC engines 
appear in the literature, diffusion and afterburning (DAB), 
and simultaneous mixing and combustion (SMC). They are 
distinguished by their thermodynamic cycle. In the DAB 
cycle, an inert rocket exhaust is used as the primary flow in 
an ejector process. Following ejector pumping, the mixed 
flow is diffused, fueled and burned subsonically in an after- 
burner. Subsequently, the gases are expanded through a con- 
vergent-divergent nozzle. This scheme results in high 
thermodynamic performance, but due to the serial nature of 
the processes, requires a long duct and suffers from low 
thrust-to-weight ratio. Any combustion during the mixing 
process reduces performance, so a stoichiometric or oxygen- 
rich rocket exhaust is generally assumed. The use of a num- 
ber of small primary nozzles can shorten the required mixing 
duct length, but the design of small, stoichiometric rocket 
nozzles raises numerous issues related to structural and ther- 
mal design and reliability. 

Some of the negative aspects of the DAB cycle are elimi- 
nated at the expense of some thermodynamic performance 
with the SMC cycle. In the SMC cycle, the rocket exhaust is 
fuel-rich and provides some fraction of the fuel required for 
combustion with the entrained airflow. The rocket and air 
streams mix and burn simultaneously, eliminating the need 
for a diffuser and afterburner duct. This process generally 
results in thermal choking where mixing is complete, fol- 


NAS A/TM— 1 999-209393 


1 



lowed by expansion through the remaining area ratio. A con- 
vergent-divergent nozzle is not required, making this cycle 
more practical for RBCC engines designed to operate also in 
scramjet mode. However, the requirement for complete mix- 
ing of the rocket and air streams still results in very long 
mixing/combustor ducts. For example, the studies of 
Dimotakis 3 indicate that, for typical conditions in an RBCC 
engine at M Q = 1 , complete mixing may be reached within an 
axial distance of between 7 and 13 duct heights (for a rectan- 
gular configuration). 

This paper presents an RBCC engine concept that is 
directed even further toward the potential for high structural 
efficiency and minimum complexity by employing a low 
speed cycle in which the rocket and ramjet streams do not 
mix. Removing the requirement for mixing can shorten the 
flowpath considerably, with a corresponding reduction in 
structural weight and wetted area. Another benefit is a reduc- 
tion in risk and complexity since a single rocket element can 
now be used without regard for mixing length. Also, higher 
thermodynamic performance is possible in other modes 
where a shorter mixing duct would reduce expansion process 
losses. 

In this new “Independent Ramjet Stream” (IRS) cycle, 
the airstream is fueled independently using the ramjet and 
scramjet mode fuel injectors located in the inlet diffuser, as 
shown in Fig. 2a. This can be accomplished upstream, since 
the air stagnation temperatures during this mode are not high 
enough to cause autoignition. The rocket will serve as a pilot 
for the fueled airstream. At the point of ignition, a flame 
propagates across the ramjet duct forming a thermal throat. 
The length of the flowpath is now determined by the flame 
propagation speed. Since the rocket is not the fuel source for 
the airstream, the rocket oxidizer-to-fuel ratio (O/F) can be 
fixed at an optimum value for best system performance. An 
additional advantage of the IRS cycle is that the fuel injec- 
tors provide the means to control the location of the thermal 
throat by adjusting the amount and radial distribution of the 
fuel injected into the airstream. 

The objective of this work is to describe a computational 
analysis tool developed for the IRS cycle, and investigate the 
performance of this propulsion mode over its operating flight 
range (0<A/ 0 <3). An additional objective is to generate 
performance maps for trajectory optimization. 

Numerical Formulation 

The ramjet and rocket streams are solved simultaneously 
using a total variation diminishing (TVD) MacCormack 

time-marching scheme 4 . It is assumed that the primary 
(rocket) and secondary streams do not mix, and therefore, 
each stream can be computed independently. Hie only inter- 
action between the primary and secondary streams occurs 
through a matched-pressure interface. 


A quasi-one dimensional approximation is used to model 
both streams. Combustion in the ramjet stream is modeled 
by a prescribed distribution of hydrogen fuel along the com- 
bustor duct as described later. Equilibrium chemistry is used 
to model the combustion process utilizing the LSENS kinet- 
ics code of Radhakrishnan 5 . The pressure matching is 
accomplished in the present work, by adding an auxiliary 
equation for the area of the primary stream, A p . 

The quasi one-dimensional Euler equations in conserva- 
tion law form for the two streams, and the auxiliary equation 
can be expressed as 


dQ.dF 
dt dx 


= S 
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The dependent variables are the mass density p, the 
momentum m (equal to p u , where u is the velocity), and 
the total energy per unit volume e . The subscripts ‘p ? and 4 s’ 
denote the primary (rocket) and secondary streams respec- 
tively. The last dependent variable is the primary stream area 
A p . The sum of the primary and secondary areas must equal 

the prescribed total engine area, A tvt , at every point. That is 


A p (x) + A s (x) = A,Jx) (3) 


The term q denotes the external rate of heat addition per 
unit mass. This term includes only the heat that is added 
externally and does not include heat released by chemical 
reactions. The variable K is a constant that controls the rate 
of change of the primary' area (0 < K < 0.5 in this study), 
and 5 is a function of the primary Mach number, M p , given 
by 


| 1 if M p > 1 

1 -1 if Mp < 1 


(4) 
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An additional equation is needed to relate the pressure to 
the flow variables p, m , and e . From the definition of static 
enthalpy per unit mass, h , we have the following relation 

h = e + p/p (5) 

where e is the internal energy per unit mass. The parameter 

(3 is introduced and defined as the ratio of static enthalpy to 
internal energy 


written as follows: 


A Q { P = 




QT = A Q/ + Q) 


a el 2) = ^[-Ae^-fVyi.-FyvA^ 1 


( 9 ) 

( 10 ) 

( 11 ) 
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q'P = a ef+e‘ n 


( 12 ) 


The total energy can be written as the sum of the internal 
energy and kinetic energy 


e/p = e + t M 2 (7) 

Combining Eqs. (5)-(7) yields the desired relationship 

p = (P-l)(e-ipH 2 ) (8) 

Equation (8) is a general expression valid for any gas mix- 
ture. The parameter P is a property of the gas and in general 
depends on its pressure, temperature and species composi- 
tion. For an ideal gas, P equals the ratio of specific heats, y. 

For a gas that is in chemical equilibrium, the parameter P is 
calculated, in this study, using the equilibrium composition 
method developed by Gordon and McBride 6 , and using the 
implementation of Radhakrishnan as given in the LSENS 
chemical kinetics code 5 . The chemical equilibrium state for 
a given initial mixture composition can be obtained when 
any two independent thermodynamic state variables are 
specified. In the present code, we compute the equilibrium 
composition based on two sets of variables: specific internal 
energy and specific volume density (UV), and temperature 
and pressure (TP). The chemical equilibrium state is com- 
puted for the TP problem by minimizing the Gibbs function, 
and for the UV problem by minimizing the Helmholtz func- 
tion. A Detailed description of the equations and solution 
procedure used to compute the chemical equilibrium state is 

given by Radhakrishnan 5 . 


Qj + ] = + + 


,< 2 ) 


1/2%+ 1/2 


(13) 


The first two steps (Eqs. 9-12) represent the standard 
MacCormack scheme where the superscripts “(1)” and “(2)” 
designate values of the function evaluated at the intermediate 

solutions Q [l) (predictor step) and Q ( ~ ] (corrector step). 
The last step (Eq. 12) represents an additional conservative 
dissipation designed in such a way that the final scheme sat- 
isfies the TVD conditions (see Ref. 7). Let 

A = (3F)/(9(?), then the vector X n J + l/2 in Eq. (13) 
denotes the matrix of eigenvectors of the flux Jacobian 
matrix A evaluated at some symmetric average of Q } and 

Qj +] denoted as l/2 . In the present work we use a 
symmetric TVD scheme for which the elements of the vector 
0" + j/ 2 , denoted by <t> 7 + 1/2 are given by: 


( Pj + 1/2 - 9 1/2 f ](«' 


+ 1/2 


0 


> + 1/2 ) 
(14) 


i _ A t i 

V 7 + 1/2 - + 1/2 


(15) 


Here a 7 + l/ 2 denotes the eigenvalues of the Jacobian 

matrix A evaluated at g - + 1/2 , and a 7 + 1/2 denotes the ele- 
ments of the vector a - + ]/2 : 

a > + !/2 = Xj+ \/2(Qj+ l ~Qj) 06) 


Numerical Method 


The function 'F is: 


The numerical method selected for solving the quasi one- 
dimensional Euler equations is the explicit, second-order, 

Total Variation Diminishing (TVD) MacCormack scheme 4 . 
The auxiliary primary area equation is integrated using a 
simple Euler method. The TVD MacCormack scheme can be 


V(z) = 


\z\ 

q 2 + e 2 ) 
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(17) 
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The entropy correction parameter, e , in Eq. (17) is taken 
to be a small positive number in the range 0.05 < e < 0.25 . A 
value of e = 0.125 was used in the present work. The lim- 
iter function 0^ + 1/2 used in this study is the following: 

e v + i/2 = minmod(a'_ ]/2> a' + 1/2 , a' +V2 ) (18) 

The expressions for the eigenvalues and eigenvectors are 
given in the Appendix. 

Heat Release Distribution 

A brief description of the combustion process in the IRS 
cycle is needed in order to understand the heat release model 
adopted in this work. In the IRS cycle, fuel injectors located 
in the inlet diffuser are used to distribute fuel into the air- 
stream. The premixed stream is then ignited by contact with 
the hot rocket plume, and a flame propagates across the sec- 
ondary stream. The combustion is completed at some point 
before the end of the ramjet duct. If the fuel is distributed in 
such a way that it does not penetrate completely across the 
duct and instead a thin layer of noncombustible gas exists 
between the rocket plume and the fueled airstream, then 
some mixing and diffusion must take place for ignition to 
occur. This radial distribution of fuel can therefore be used to 
control the point of ignition, and therefore the thermal throat 
location. 

In the present quasi-one dimensional analysis, we mode] 
the combustion process by specifying an ignition point, x s q , 

an end of heat release point, x™ , and a distribution function 
between them. For ideal gas calculations, the heat release 
distribution is specified by the source term q{x) , where the 

total heat added per unit time is given by jp s qdV . For equi- 
librium flow calculations, the heat release distribution is 
specified as follows. Consider a case where the fuel is hydro- 
gen. Then, we define a new chemical species named “inert 
hydrogen” which has identical properties as “real” hydrogen 
with the exception that inert hydrogen does not participate in 
any chemical reactions. The total hydrogen concentration at 
any given point is then given by the sum of inert and real 
hydrogen. By prescribing the axial distribution of real hydro- 
gen along the flow domain we can specify the heat release 
distribution, and at the same time introduce a combustion 
efficiency (by specifying a nonzero value for the inert hydro- 
gen at the end of the combustion process). For the RBCC 
computations, the following cubic function was used to 
express the real hydrogen mole concentration 



where r\ c is the combustion efficiency, <|> is the equiva- 
lence ratio, = x e q - x*j is the length of the heat release 
zone, and z is the nondimensional distance, 
z = (x-x^)/(Ax q ). 

Results 

The accuracy of the method is first demonstrated by pre- 
senting computations for two benchmark cases, involving 
reacting and nonreacting flows. They are presented below. 

Benchmark test cases 

The first case considers the nonreacting flow of an ideal 
gas in a divergent nozzle. We use the unsteady quasi-one 
dimensional Euler equations to obtain the steady-state solu- 
tion for an inflow Mach number M l = 1.25 and a back 

pressure p b = 1, 87466/? x , where p ] is the static pressure at 

the inflow plane. The area distribution A(x) of the divergent 
nozzle considered is given in Fig. 3, which compares the 
CFD result with the exact analytical solution. For these flow 
conditions, a normal shock forms at the midpoint of the noz- 
zle (jc = 0.5). The CFD solution was computed on a 200 
point grid and demonstrates the accuracy of the TVD 
scheme, which leads to sharp stationary shock transitions, 
over two grid points. 

The second test case considers thermal choking in a clas- 
sical Rayleigh flow, which consists of frictionless, one 
dimensional flow in a constant area duct with heat addition. 
Calculations were performed for an inflow gas mixture of 
hydrogen-air having a fuel-air ratio of f /a = 0.00152 , and 
a sufficiently low exit pressure to cause thermal choking. If 
we assume that the hydrogen-air mixture behaves as an ideal 
gas, then an exact analytical solution exists for this problem 
(Ref. 9). The total temperature ratio T t2 /T tl between the 
exit and entrance to the duct can then be easily computed to 
be T tl /T tX = 1.4464 for this fuel-air ratio. For a choked 

flow (A/ 2 = 1.0), the inflow Mach number, M { , is analyt- 
ically computed to be M x = 0.5 , and the corresponding 
total pressure loss is p t2 ^P t \ = 0.8976 . 

Figure 4 shows the CFD results for the above conditions 
assuming an ideal gas. The external heat is added from sta- 
tions x - 0.25 to x = 0.75 . The computed inflow Mach 
number is M j = 0.5004 , and the computed total tempera- 
ture and total pressure ratios are T tl /T tX = 1.4452, and 
p t2 /p t \ = 0.8978 respectively, which are in excellent 
agreement with the analytical solution. 
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The calculation was repeated for the same fuel-air ratio 
but without the assumption of ideal gas, and considering 
instead a real gas in chemical equilibrium. The heat release 

zone considered was from x s q * = 0.25 to x q = 0.75 , using 

a cubic function slightly different from the one given in Eq. 
19, and different also from the distribution assumed in the 
ideal gas calculation. (We should point out that for this case, 
the inflow and outflow states are independent of the assumed 
heat release distribution.) The results are shown in Fig. 5, 
which compares the pressure, temperature and Mach number 
distribution for equilibrium and ideal gas flows. The com- 
puted inflow Mach number for the equilibrium calculation is 
slightly higher than the theoretical value of = 0.5 , 
probably due to a small amount of dissociation at the exit 
plane which results in a slightly lower amount of effective 
heat release. The assumption of a constant specific heat in 
the ideal gas calculation may contribute also to the small dif- 
ferences observed between the ideal gas and equilibrium 
flow calculations. 

Analysis of a Trailblazer RBCC engine at low speed 

We investigated the performance characteristics of the 
IRS cycle for a specific Trailblazer configuration over the 
flight Mach number range of 0 < A/ 0 < 3 . All computations 
were carried out with equilibrium chemistry. 

Referring to the control volume shown in Fig. 2b, the net 
thrust is calculated by 


j — 1 neI 

SP ~ *p + (”>w).s 


( 24 ) 


Here, w p is the propellant weight flow rate for the pri- 
mary stream, and (wj uel ) s is the fuel weight flow rate for the 
secondary stream. 


The rocket element is designed to operate at a constant 
O/F ratio but at variable chamber pressure. Therefore, we 
investigated the performance of the IRS cycle at different 
chamber pressures, p c , in the range 100 < /? c .//? 0 < 350 , 

where p 0 is the free- stream pressure. The design parameters 
for the configuration studied are given below: 


Capture area (per pod) A c = 18094.7 in 2 
Ramjet duct inlet total area ratio A$/ A c = 0.4 
Ramjet duct exit area ratio A 6 /A c = 1 .25 

Ramjet duct nondimensional length L/A c = 2.1452 

Ramjet duct inlet secondary area ratio A 3 /A c =0.3318 

Rocket O/F ratio = 6.0 

Rocket throat area ratio A /A c = 0.0068 

* 

Rocket exit area ratio A$ p /A - 10.0 


F — F — F F 

* net r gross r ram r inlet 

(20) 

The first term on the right hand side of Eq. (20) is the 
gross thrust, and is computed as follows 

F ideal - P9 W 9^9 

(21) 

Z7 F F 

r gross ^ fg r ideal 

(22) 

where C. is the nozzle gross thrust coefficient. The sec- 

ond term on the right hand side of Eq. (20) is the 
and can be written as 

ram drag. 

F ram = P\ u U\ 

(23) 


The final term in Eq. (20) represents the inlet drag which 
includes the contributions from spillage, cowl lip and 
diverter drag, and is obtained from the two-dimensional CFD 

calculations of the inlet reported by DeBonis et a/ 8 . 


Maximum combustion efficiency r| r = 0.95 

The calculations assumed that the rocket (primary) flow 
was frozen from the throat downstream, and that the second- 
ary flow is in chemical equilibrium. The heat release distri- 
bution was based on the cubic function given in Eq. (19). 
The CFD computations were carried out from station 3 to 
station 6 (see Fig. 2b). An isentropic expansion to free- 
stream conditions was assumed from station 6 to station 9. 
Expansion losses were included in the nozzle gross thrust 
coefficient Cj g . The values assumed for Cj g are given in 
Table 1. 

The boundary conditions were specified as follows. At 
the inflow plane (station 3 in Fig. 2b), the primary stream 
pressure, temperature, Mach number and gas properties were 
specified. For the secondary stream, the total pressure, total 
temperature and fuel-air ratio was specified. Total pressure 
values were based on the inlet recovery of Ref. 8. At the out- 
flow boundary (station 6), the outflow pressure was imposed 
if the flow was subsonic, and all the variables were extrapo- 
lated if the flow was supersonic. 


The net specific impulse, I sp , is then calculated by 


Since the maximum amount of air flowing into the 
engine is known at every Mach number from the CFD inlet 
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calculations 8 , the fuel-air ratio in the secondary stream, and/ 
or the heat release distribution were adjusted to match this 
given airflow. 

Figures 6 and 7 show typical flow properties distribution 
for the primary and secondary streams along the ramjet duct. 
The distance x has been normalized by the ramjet duct 
length, thus x = 0 corresponds to station 3, and x = 1 cor- 
responds to station 6. This figure shows that the primary and 
secondary pressures are matched within a couple of grid- 
points. At the Mach-1 condition, the rocket flow is underex- 
panded and therefore must expand rapidly to match the 
secondary stream pressure. As a result, the primary Mach 
number increases and the static temperature decreases. Fur- 
ther downstream, the primary Mach number continues to 
increase gradually through the expanding ram jet duct at the 
same time that the temperature decreases. At the Mach-2 
condition, the rocket is overexpanded and its pressure must 
be brought up to that of the secondary stream. The rapid rise 
in primary pressure and temperature (and a decrease in Mach 
number) is the one dimensional approximation of an oblique 
shock. The pressure in the ramjet duct remains high until the 
heat release causes the pressure to fall. For both of these 
cases, the secondary stream thermally chokes at a location 
corresponding approximately to the end of the heat release 
(maximum secondary flow temperature). 

Figure 8 shows a typical convergence history for the pri- 
mary and secondary flows. For all cases presented, the calcu- 
lations were stopped when both norms were smaller than 

10' 6 . Experience showed that continuing the computations 
far beyond this convergence level resulted in only negligible 
changes in the performance parameters (less than half a per- 
cent). 

Figure 9 shows the net specific impulse variation with 
free-stream Mach number for various rocket chamber pres- 
sure ratios. The effect of ramjet cycle efficiency is evident 
above Mach 1.5. Also apparent is the trend toward higher 
specific impulse as the rocket chamber pressure is reduced. 
The dip in specific impulse at Mach 1 is due to inlet drag. 

Figure 10 shows the variation of net thrust coefficient, 
defined as C { = F net /(p 0 A c ) , with Mach number. Above 
Mach 1.5, the percentage of thrust produced by the ramjet 
increases rapidly. Therefore, as the vehicle accelerates, the 
rocket can be throttled down to increase the efficiency. When 
the ramjet thrust alone is sufficient for adequate acceleration, 
the rocket can be completely shut off. Therefore, a smooth 
transition from rocket propulsion at lift-off to ramjet propul- 
sion at some design supersonic Mach number can be accom- 
plished. 

Figure 1 i shows the ratio of the secondary to primary 
flow, P , defined as 


0 = (25) 

For clarity, the plot is presented only for the highest and 
lowest chamber pressures considered in the study. The 
results indicate that at high supersonic speeds, and low 
rocket chamber pressure, the secondary flow can be up to 25 
times larger than the rocket flow. 

At each Mach number between Mach 1 and Mach 2.5, a 
different fuel-air ratio was needed for each value of the 
rocket chamber pressure in order to keep the airflow fixed at 
the inlet critical value. More fuel is required at the lower 
chamber pressures. Figure 12 shows the required equiva- 
lence ratio as a function of Mach number for the various 
chamber pressures. Below a Mach number of 1.0, the air- 
stream is not fueled since there is not enough ram pressure 
for efficient ramjet operation. At Mach 2.5 and above, the 
fuel-air ratio reaches stoichiometric proportions for all 
chamber pressures. From Mach 1 to Mach 2, the same heat 

distribution given in Eq. 19 was used with x q = 0.2, and 
x q = 0.95 The thermal throat for all these cases was estab- 
lished just ahead of x ™ , going from x* = 0.94 at Mach-1 

to x* = 0.90 at Mach-2. For Mach numbers above 2.5, 
where the secondary stream was already at stoichiometric 
proportions, it was necessary to move the thermal throat 
upstream in order to match the prescribed airflow. This was 
accomplished by moving the ignition point upstream and 
shortening the heat release zone. The thermal throat location 
was approximately x* = 0.77 for Mach 2.5, and x* = 0.45 
for Mach 3. 

Figures 13 and 14 show the variation of specific impulse 
and net thrust coefficient with equivalence ratio at various 
Mach numbers, and for a rocket chamber pressure of 
p c /po = 200 . The inlet drag was assumed to be constant at 
each Mach number. The figures show that adding fuel at a 
Mach number of 0.5 does not generate any significant addi- 
tional thrust. The benefits of adding fuel to the secondary 
stream begin to become evident at Mach 1 . At higher Mach 
numbers, the specific impulse variation with equivalence 
ratio shows a rapid initial increase followed by a more grad- 
ual increase at higher values of the equivalence ratio. From 
an overall vehicle performance point of view this suggests 
that the optimum equivalence ratio should be near the “knee” 
of the curve (<j> ~ 0.6) since adding additional fuel results in 
only small increases in performance at the cost of larger, and 
heavier fuel tanks. 

Figure 15 compares the specific impulse for the IRS and 
SMC cycles from lift-off to Mach 3, along a typical Trail- 
blazer flight trajectory. The SMC calculations were carried 
out using RAMSCRAM 1 *, a control-volume type cycle anal- 
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ysis code for ramjet and scramjet engines, including ejector 
modeling. The SMC calculations were carried out using the 
same rocket chamber pressures (listed in Fig. 16), and the 
same values for the various efficiency parameters (i.e., inlet 
recovery, combustion efficiency, nozzle gross thrust coeffi- 
cient, etc.). In addition, the SMC calculations assumed com- 
plete mixing and chemical equilibrium at the thermal throat, 
and assumed no inlet losses at sea level static conditions. 
Similar to the calculations for the IRS cycle, the SMC analy- 
sis assumed flow expansion to atmospheric pressure at sta- 
tion 9 (see Fig. 2b). 

The IRS cycle is slightly less efficient than its counter- 
part SMC cycle up to around Mach 3, where the two curves 
cross over. The differences in net specific impulse are in gen- 
eral less than 10%, with a maximum of 17% occurring at 
Mach 1. This difference in low speed performance would 
result in a 2-3% change in the launch vehicle final mass. 
Advantages in weight and simplicity of the IRS cycle may 
more than compensate for the lower cycle efficiency. 

Conclusions 

The numerical model developed in this study is an effi- 
cient approach for analyzing the performance of the IRS 
cycle. Computations of benchmark test cases for reacting 
and nonreacting flows, carried out with the present numerical 
method, were in excellent agreement with exact analytical 
results. Shock waves were sharply captured over two grid 
points. 

The results obtained for the IRS cycle show that at sub- 
sonic speeds, there is no benefit in adding fuel to the air- 
stream. The airflow, however, prevents the rocket plume 
from over-expanding. Above Mach 1, the airflow can be 
fueled and burned to generate ramjet thrust. The fraction of 
the thrust generated by the ramjet part of the combined cycle 
increases rapidly with Mach number, and the rocket can be 
progressively throttled down to obtain a higher specific 
impulse. When the ramjet thrust alone is sufficient, the 
rocket can be completely shut off. Therefore, a smooth tran- 
sition from all-rocket at lift-off to pure ramjet at some design 
supersonic Mach number is naturally attained. 

An optimum value of the equivalence ratio can be 
defined at each Mach number, above which only small 
increases in thrust and specific impulse occur. Operating 
above the optimum equivalence ratio may not be justified 
due to hydrogen tank size considerations. 

Low speed operation of an RBCC engine requires the 
ability to control the inlet diffuser exit pressure in order to 
keep the inlet operating at optimum conditions through the 


entire Mach number range. The IRS cycle can achieve this 
by adjusting the amount and radial distribution of fuel 
injected into the secondary stream, thereby controlling the 
location of the thermal throat. 

Comparisons between the unmixed IRS cycle and the 
fully mixed SMC cycle showed that the IRS cycle perfor- 
mance was slightly lower than that of the SMC cycle. How- 
ever, when other advantages of the IRS cycle, such as high 
structural efficiency, minimum system complexity and con- 
stant rocket O/F ratio are considered, the IRS cycle may be 
the best choice for an RBCC single-stage-to-orbit vehicle. 
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Appendix 


1 

^(6, + u/a) 

~(ub 2 + 1/a) 

k 

jt' = 

l-*i 

b 2 u 

-b 2 


-u/a) 

-^(ub 2 - 1/a) 

k 

where b l = p p /a 2 

, and b 2 = p e /a 2 



For the quasi-one dimensional Euler equations, the flux 
Jacobian matrix A is given by 


A = 


0 1 0 
P p -U u(2-p e ) p e 
u(p p -H) H-u 2 p e u(l + p e ) 


(26) 


Here H is the total enthalpy per unit mass 

H = (e + p )/ p (27) 

Note that the relation between pressure and the flow 
variables involves P which itself is a function of density and 
energy, and therefore complicates the computation of the 
partial derivatives in the Jacobian matrix. In the present work 
we utilize expressions for these partial derivatives obtained 
for the more general case of finite rate chemistry. The result- 
ing expressions are (see Ref. 10): 

Pe = J (28) 


Pp = RT(\+p e ) + p e (u 2 -H) (29) 


Here R and c v denote the specific gas constant, and spe- 
cific heat at constant volume for the mixture at the equilib- 
rium conditions. The eigenvalues and eigenvectors of the 
Jacobian matrix are: 


A = Diagfw-tf, «, u + a] (30) 

where a is the frozen speed of sound defined by 

(31) 


a = Pp + Pe( H -u 2 ) 


X = 


1 


u- a 


1 

u + a 


(H-ua) (u -p p /p e ) (H + uc) 


( 32 ) 


Table 1: Nozzle gross thrust coefficient. 


M o 

Qg 

0.0 

0.95 

0.5 
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0.90 
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0.925 
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0.95 
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0.95 

3.0 

0.95 


« 


NASA/TM— 1999-209393 


8 





Figure 1 . Trailblazer SSTO Vehicle. 
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Figure 2. Cross-sectional view of the Trailblazer engine.; (a) schematic of the Independent Ramjet 
Stream (IRS) propulsion mode, (b) Control Volume around the RBCC engine. 
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Figure 3. Numerical and exact solutions of ideal gas flow in a 
divergent nozzle. 



Figure 4. Variation of Mach number, total pressure and total 
temperature in Rayleigh Flow (ideal gas). 



Figure 5. Comparison between equilibrium and ideal gas cal- 
culations for Rayleigh Flow. (PI = 1.0 atm., T 1=404.7 K are 
reference values.) 


Figure 6. Pressure, temperature and Mach number distribu- 
tion for A/ 0 = 1 , p c / p 0 = 200 . Subscripts are: p- primary, 
s- secondary, 0- free- stream. 



X 


Figure 7. Pressure, temperature and Mach number distribu- 
tion for M 0 = 2, p c /p 0 = 200. 
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Figure 8. Convergence history for primary and secondary 
flows. 
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Figure 9. Net specific impulse as a function of free stream Fjgure , 2 Equivalence ratio as a function of free stre am 

Mach number for various rocket chamber pressure ratios. Mach number for various rocket chamber pressure ratios 




Figure 10. Net thrust coefficient as a function of free stream Figure 13. Net specific impulse as a function of equivalence 

Mach number for various rocket chamber pressure ratios. ratio for various free-stream Mach numbers p c /po = 200. 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Free-stream Mach number, M 0 

Figure 1 1. Ratio of secondary to primary flow as a function 
of free stream Mach number for two rocket chamber pres- 
sure ratios. 



Figure 14. Net thrust coefficient as a function of equivalence 
ratio for various free-stream Mach numbers, p c /p 0 = 200. 
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Free-stream Mach number, M 0 


Figure 15. Net specific Impulse comparison for IRS and 
SMC cycles along a typical Trailblazer trajectory from lift- 
off to Mach 3. 
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